function specifications = pick_var_fn(model, settings, spec_id)
%Modified by Viet Hoang Dinh

% Function for randomly choosing DGPs from the encompassing model
    % Warning: a "specification" here is a "DGP" in the paper

    % prepare
    
    rng(spec_id, 'twister');

    specifications = settings.specifications;
    
    n_y = model.n_y;
    
    manual_var_select = specifications.manual_var_select;
    random_select = specifications.random_select;
    random_category_range = specifications.random_category_range;

    % randomly choose
    
    if random_select == 0

        var_select = manual_var_select;

    else
        
        random_n_spec         = specifications.random_n_spec;
        random_n_var          = specifications.random_n_var;
        random_fixed_var      = specifications.random_fixed_var;
        random_fixed_pos      = specifications.random_fixed_pos;
        random_fixed_var2     = specifications.random_fixed_var2; %Position of interested variable e.g. IP or Unem

        random_category_setup = specifications.random_category_setup;

        var_select = nan(random_n_spec, random_n_var);
        var_select(:, 1) = random_fixed_var;
        fixed_category = sum(random_fixed_var > random_category_range(:,2)) + 1; % category of fixed variable
        
        for i_spec = 1:random_n_spec % draw for each specification
            
            for iv = 1:(random_n_var - 1) % draw non-fixed variables from other categories
                
                % check if draw from certain categories
                if iv > length(random_category_setup)
                    feasible_range = 1:n_y;
                else
                    feasible_range = [];
                    for icat = random_category_setup{iv}
                        feasible_range = [feasible_range, random_category_range(icat,1):random_category_range(icat,2)];
                    end
                end

                % draw this random variable
                while true
                    
                    drawn_pos = randi(length(feasible_range));
                    drawn_ID = feasible_range(drawn_pos);
                    
                    if all(var_select(i_spec, 1:iv)~=drawn_ID) % add this new variable only if no same variable exists
                        var_select(i_spec, iv + 1) = drawn_ID;
                        break;
                    end
                   
                end

            end
            
            % randomly permute non-fixed variables
            var_select(i_spec, 2:random_n_var) = var_select(i_spec, 1 + randperm(random_n_var - 1));
            
        end
        
        
        % put fixed variable at the fixed position
        var_select(:, [1 random_fixed_pos]) = var_select(:, [random_fixed_pos 1]); % put fixed var at fixed position
        
    end

%New code: rearrange the order of variables
var_select2=nan(random_n_spec, random_n_var);
for i_spec = 1:random_n_spec
    var_select_row=var_select(i_spec,:);
    var_select_row2=var_select_row;
    for i_var=1:n_y
        if var_select_row(:,i_var)==random_fixed_var2       
            var_select_row2(:,i_var)=[];
            var_select2(i_spec,:)=[var_select_row2(:,1) random_fixed_var2 var_select_row2(:,2:end)];            
        end
    end
    
end

   var_select=var_select2;
   n_spec = size(var_select, 1);
   n_var  = size(var_select, 2);

    
    specifications.var_select = var_select;
    specifications.n_spec     = n_spec;
    specifications.n_var      = n_var;


end